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■ ABSTRACT 

^ \ X-ray photons emitted from the surface or atmosphere of a magnetized neu- 

0^ ■ tron star are highly polarized. However, the observed polarization may be 

o ■ 

■ modified due to photon propagation through the star's magnetosphere. For 
CO ' photon frequencies much larger than the typical radio frequency, vacuum bire- 

o ; 

^\ ^ fringence due to strong-field quantum electrodynamics dominates over the 

, plasma effect. We study the evolution of photon polarization in the magne- 

> • 

•'^ ■ tized QED vacuum of a neutron star magnetosphere, paying particular atten- 

\^ ' tion to the propagation effect across the quasi-tangential (QT) point, where 

the photon momentum is nearly aligned with the magnetic field. In agreement 
with previous studies, we find that in most regions of the magnetosphere, the 
photon polarization modes are decoupled due to vacuum birefringence, and 
therefore a large net linear polarization can be expected when the radiation 
escapes the magnetosphere. However, we show that the X-ray polarization 
may change significantly when the photon passes through the QT region. 
When averaging over a finite emission area, the net effect of QT propagation 
is to reduce the degree of linear polarization; the reduction factor depends on 
the photon energy, magnetic field strength, geometry, rotation phase and the 
emission area, and can be more than a factor of two. We derive the general 
conditions under which the QT propagation effect is important, and provide 
an easy-to-use prescription to account for the QT effect for most practical 
calculations of X-ray polarization signals from magnetic neutron stars. For a 
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neutron star with a dipolc magnetic field, the QT effect can be important for 
emission from the polar cap for certain magnetic field and energy ranges, and 
is negligible for emission from the entire stellar surface. 

Key words: plasmas - polarization - waves - star: magnetic fields - pulsars: 
general - X-rays: stars 

1 INTRODUCTION 

Thermal, surface emission from neutron stars (NSs) has the potential of providing invaluable 
information on the physical properties and evolution of NSs (equation of state at super- 
nuclear densities, cooling history, magnetic field, surface composition, different populations; 
see, e.g., Yakovlev & Pethick 2004; Harding & Lai 2006). With X-ray telescope such as Chan- 
dra and XMM-Newton, the last decade has seen significant observational progress, revealing 
the surface magnetic field geometry of isolated pulsars with phasc-resolvcd spectroscopy, and 
constraining the cooling physics from thermal emission of young NSs in supernova remnants 
(e.g., Kaspi et al. 2006). In addition, thermal emission from seven isolated, radio-quiet NSs 
has been studied in detail, revealing absorption features in their spectra in many cases (see, 
e.g., van Kerkwijk & Kaplan 2007; Kaplan 2008). 

It has been recognized that in addition to imaging, timing and spectroscopy. X-ray po- 
larimetry provides a new way to study many high-energy astrophysical sources, particularly 
magnetic NSs. Recent advances in detector technology suggest that polarimetry study of 
X-ray sources holds great promise in the future (e.g., Costa et al. 2008; Swank et al. 2008). 

The surface emission from magnetized NSs (with B ^ 10^^ G) is highly polarized (e.g., 
Gnedin & Sunyaev 1974; Meszaros et al. 1988; Pavlov & Zavhn 2000) for the following reason. 
In the magnetized plasma that characterizes NS atmospheres. X-ray photons propagate in 
two normal modes: the ordinary mode (0-mode, or ||-mode) is mostly polarized parallel 
to the k-B plane, while the extraordinary mode (X-mode, or _L-mode) is mostly polarized 
perpendicular to the k-B plane, where k is the photon wave vector and B is the external 
magnetic field (e.g., Meszaros 1992). This description of normal modes apphes under typical 
conditions, when the photon energy E is much less than the electron cyclotron energy 
Ebc = heB/{mec) = 11.6Si2keV [where S12 = ^/(lO^^G)], E is not too close to the 
ion cyclotron energy Esi = 6.3 Bi2{Z/A) eV (where Z and A are the charge number and 
mass number of the ion), the plasma density is not too close to the vacuum resonance (see 
below) and 9b (the angle between k and B) is not close to zero. Under these conditions, the 
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X-mode opacity (due to scattering and absorption) is greatly suppressed compared to the 
0-mode opacity, kx ~ {E/EbsYi^o (e.g. Lodenquai et al. 1974; Meszaros 1992; Potekhin & 
Chabrier 2003). As a result, the X-mode photons escape from deeper, hotter layers of the 
NS atmosphere than the 0-mode photons, and the emergent radiation is linearly polarized 
to a high degree (e.g., Pavlov & Zavhn 2000; Ho & Lai 2001,2003; van Adelsberg & Lai 
2006). Measurements of X-ray polarization, particularly when phase- resolved and measured 
in different energy bands, could provide unique constraints on the NS magnetic field strength 
and geometry. 

It has long been predicted from quantum electrodynamics (QED) that in a strong mag- 
netic field the vacuum becomes birefringent (e.g., Schwinger 1951; Adler 1971). While this 
vacuum polarization effect makes the photon index of refraction deviate from unity only when 
B ^ 300Sq, where Bq = mlc^/{eh) = AAU x 10^^ G is the critical QED field strength, it can 
significantly affect the spectra of polarization signals from magnetic NSs in more subtle way, 
at much lower field strengths (see section 2 of Lai & Ho 2003a for a qualitative explanation). 
In particular, the combined effects of vacuum polarization and magnetized plasma gives rise 
to a "vacuum resonance", at which the contributions from these two effects (plasma and 
vacuum polarization) to the dielectric tensor "compensate" each other (Gnedin et al. 1978; 
Meszaros & Ventura 1979; Pavlov & Gnedin 1984; Lai & Ho 2002). A photon may convert 
from the high-opacity mode to the low-opacity one and vice verse when it crosses the vacuum 
resonance region in the inhomogcncous NS atmosphere. For B > 7 x 10^"^ G, this vacuum 
resonance phenomenon tends to soften the hard spectral tail due to the non-greyness of the 
atmospheric opacities and suppress the width of absorption lines, while for S ^ 7 x 10^^ G, 
the spectrum is unaffected (see Lai & Ho 2002,2003a; Ho & Lai 2003; van Adelsberg & Lai 
2006). 

The QED-induced vacuum birefringence influences the X-ray polarization signals from 
magnetic NSs in two different ways, (i) Photon mode conversion in the NS atmosphere: Since 
the mode conversion depends on photon energy and magnetic field strength, this vacuum 
resonance effect gives rise to a unique energy-dependent polarization signal in X-rays: For 
"normal" field strengths {B ^ 7 x 10^^ G), the plane of linear polarization at the photon 
energy E <> IkeV is perpendicular to that at ^ 4keV, while for "superstrong" field 
strengths {B > 7 x 10^"^ G), the polarization planes at different energies coincide (Lai & Ho 
2003b; van Adelsberg & Lai 2006). (ii) Polarization mode decoupling in the magneto sphere: 
The birefringence of the magnetized QED vacuum decouples the photon polarization modes. 
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so that as a polarized photon leaves the NS surface and propagates through the magneto- 
sphere, its polarization direction follows the direction of the magnetic field up to a large 
radius (the so-called polarization limiting radius). The result is that although the magnetic 
field orientations over the NS surface may vary widely, the polarization directions of the 
photon originating from different surface regions tend to align, giving rise to large observed 
polarization signals (Heyl & Shaviv 2002; van Adelsberg & Lai 2006). 

In this paper we examine in detail the photon polarization evolution in the magnetized 
QED vacuum of a NS magnetosphere [Effect (ii) in the last paragraph] . We are particularly 
interested in photon propagation across the Quasi-tangential Region (QT region): As the 
photon travels through the magnetosphere, it may cross the region where its wave vector 
is aligned or nearly aligned with the magnetic field (i.e., is zero or small). In such a 
QT region, the two photon modes (|| and ± modes) become (nearly) identical, and mode 
coupling may occur, thereby affecting the polarization alignment. Our previous analytical 
(or semi-analytical) treatment of the alignment effect (Lai & Ho 2003b; van Adelsberg & 
Lai 2006) focused on the region far away from the NS surface, thus did not not include the 
QT region (which typically lies within a few stellar radii). The numerical ray integration 
presented in Heyl & Shaviv (2002; see also Heyl et al. 2003) should in principle have included 
such QT region, but no systematic characterization of the QT propagation effect on the final 
photon polarization was attempted there. As we show in this paper, polarization evolution 
through the QT region is sufficiently subtle (e.g. the effect varies on small length scales 
across the emission region) that a careful examination of its effect is necessary. The purpose 
of this paper is to study the evolution of high-energy (X-ray) photon polarization in NS 
magnetospheres and to quantitatively assess the QT propagation effect. 

The remainder of our paper is organized as follows: Section 2 summarizes the basic 
equations for studying photon polarization evolution in magnetized QED vacuum. In section 
3 we examine the general behavior of the polarization evolution across a QT region in a 
generic magnetic geometry. In section 4 we present detailed calculations in the case of dipole 
magnetic field and consider emissions from both the polar cap and the other regions of the 
NS surface. We provide a simple prescription (see section 4.3) for including the QT effect in 
the calculations of the observed polarization signals. In section 5 we discuss the implications 
of our results for the X-ray polarization signals from magnetic NSs and the prospect of using 
X-ray polarimetry to probe strong-field QED. 
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2 POLARIZATION EVOLUTION IN HIGHLY MAGNETIZED QED 
VACUUM: EQUATIONS 

The magnetospheres of pulsars and magnetars consist of relativistic electron-positron pairs 
streaming along magnetic field lines. The Lorentz factor 7 of the streaming motion and the 
plasma density N are uncertain. For the open field line region of radio pulsars, pair cascade 
simulations generally give 7 ~ 10^ — 10^ and rj = N/Nqj ~ 10^ — 10^ (e.g., Daugherty 
& Harding 1982; Hibschman & Arons 2001; Medin & Lai 2009), while recent theoretical 
works suggest that the corona of magnetars consist of pair plasma with 7 up to 10^ and 
77 ~ 2 X 10^ (i?*/r) (where i?* is the stellar radius; Thompson et al. 2002; Beloborodov & 
Thompson 2006), where A^gj = {^B)/{2Trec) is the Goldreich- Julian density. In general, 
both the plasma and vacuum polarization affect the photon modes in the magnetosphere. 
For a given photon energy E, 7 and B, the vacuum resonance occurs at the density (Wang 
& Lai 2007) 

AV = 5.80 X 10^^ Bl^El-fl (1 - (3 cosesf F{b) cm'^ (2.1) 

where p = -yT^^T/Y, 9b is k-B angle, = £;/(lkeV), S13 = B/{1Q^^G), 73 = 7/10^ 
F{b) is equal to unity for b — B/Bq <^ 1 and is at most of order a few for B ^ 10^^ G. For 
the typical photon energy of interest in this paper, the magnetosphere plasma density N is 
much less than Ny, and vacuum birefringence dominates over the plasma effect. To put it 

in another way, at given B, 7 and density N (or rj), we can define the vacuum resonance 
photon energy: 

Ey = 3.48 X 10"^° [Pr^^rs^Ts'^ (1 - /9cos^b)~^ F-'^f^^ keV. (2.2) 

where Pi is the NS spin period in units of 1 second. Throughout this paper, wc shall be 
interested in photon energies E ^ Ey, so that the wave modes are determined by the 
vacuum polarization effect. 

In this section we summarize the key equations for studying the polarization evolution 
of X-rays in NS magnetospheres. 

2.1 Wave Modes 

The dielectric tensor and the inverse permeability tensor of a magnetized QED vacuum take 
the form: 

e^al + qBB, /j,'^ ^ al + mBB, (2.3) 
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where / is the unit tensor and B = B/B is the unit vector along B. In the low frequency 
limit, E — hcu <^ meC^, the general expressions for the vacuum polarization coefficients a, q, 
and m are given in Adler (1971) and Heyl & Hernquist (1997). For S < Sq = mlc^/{eh) = 
4.414 X 10^^ G, they are given by 

a = l-25v, ? = 75v, m = -4(5v, (2.4) 
with 



Here a-p = e^/hc = 1/137 is the fine structure constant. For B ^ Bq, simple expressions 
for a, q, m are given in Ho & Lai (2003) (see also Potekhin et al. 2004 for general fitting 
formulae) . 

Using the relations D — e -E, B — jx H and the Maxwell equations, we obtain the 
equation for plane waves with E oc e*('°"^~'^*) 
1 



-eij + 



kikj - Sij {k X B)i{k x B)j \ Ej - 0, (2.6) 



where n = ck/ui is the refractive index and k — k/k. The dielectric tensor are given by 
eq. (2.3). In the coordinate system xyz (with k along the 2;-axis and B in the x-z plane, 
such that k X B = —smOsy), we can solve equation (2.6) to obtain the two eigenmodes: 
the ||-mode (or ordinary mode, polarized in the k-B plane) and ±-modes (or extraordinary 
mode, polarized perpendicular to the k-B plane). The refractive indices and polarization 
states of these two modes are 

/ (; \ V2 

— \Ey/E^\^ ^oo. (2.7) 

The difference between the two refractive indices is (for q, \m\ <^ 1) 

An = — nj_ ~ -(g + m) sin^ 6b- (2.8) 



2.2 Mode Evolution Equations and the Adiabatic Condition 

A general polarized electromagnetic wave with frequency cu traveling in the z-direction can 
be written as a superposition of the two modes: 

E{z)^Ai\{z)Eii{z) + A4z)E4z), (2.9) 
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Note that both Ay, A± and E^y, E± depend on z. Substituting equation (2.9) into the wave 
equation 



2 

V X (^Ij,-^ -V X E) ^ —e- E, (2.10) 

we obtain the amphtude evolution equations (see Adelsberg & Lai 2006) 
/ 



. d 



\ 




(2.11) 



Ay \l -{u/c)An/2 i^'^ 
A± j y -lip's (u;/c)An/2 

where ip'^ — dc/^fi/dz. Here </?b is the azimuthal angle of B in the fixed frame XYZ with 
the Z-axis parallel to the line of sight k (Note that the xyz-frame introduced in section 2.1 
rotates around the Z-axis since B changes along the ray). In deriving equation (2.11), we 
have used the geometric optical approximation \dA/dz\ <C k\A\. 
The condition for the adiabatic evolution of wave modes is 

^ » 1. (2.12) 

Here Fad is the adiabaticity parameter which changes along the photon ray. In the adiabatic 
limit, the photon modes are decoupled, and the photon will always stay in the initial mode 
although the mode itself changes along the ray following the direction of B^^ (the projection 
of B in the XY plane). Substituting eq. (2.8) into eq. (2.12), we have 



ad 



ad 



where 



(g + m)a;sin 6b 



Ac If'B 



~ 1.0 X 10^ EiBls^^^^F{b), (2.13) 



is equal to unity for b = B/ Bq ^ 1 and is at most of order a few for B < 10^^ G (sec Fig. 1 
of Ho & Lai 2003) — wc will use F ~ 1 hereafter. The unit of (p'^ is km~^ and the factor 
sin^ 9b/<p'b is determined by magnetic field geometry along the ray. 

Once the mode amplitudes A\\ and A± are known, we can obtain the Stokes parameters 
in the fixed XYZ-ir&me from (see van Adelsberg & Lai 2006) 

Q = cos2(^B(|Ay|2 - \A±\'^) - 2sin2(^B3?e(AyAl), 

U = sm2ipBi\A\\\^ -\A±\'^) + 2cos2ipB^e{AiiAl), 

V = 2^m{A\\Al). (2.15) 
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Alternatively, we can use equation (2.11) to obtain the evolution equation for the Stokes 
parameters (see van Adelsberg & Lai 2006): 

dQ/dZ = {(J / c)AnV sin 2ipB, 
dU/dZ = -{iu/c)AnV cos2(pB, 

dV/dZ = -(u;/c)An(Qsin2^s - C/cos2(^s). (2.16) 
and dl/dZ = 0. 



3 POLARIZATION EVOLUTION ACROSS THE QUASI-TANGENTIAL 
REGION 

Here we examine the general behavior of the polarization evolution when a photon crosses 
the quasi-tangential (QT) region. An exact tangential point is where the photon wave vector 
k is aligned with B, or 9b — 0. Across the tangential point, the azimuth of the magnetic 
field, (fB, changes by 180°. For a general magnetic field geometry, not all photon rays have 
such an exact tangential point. However, there exits a quasi-tangential point where the k-B 
angle 6b reaches a local minimum. Without loss of generality, the magnetic field around the 
QT point can be modeled as 

Bx = ^s, BY = eB, (3.17) 

in the fixed XYZ-bame with Z \\ k. Here B is the magnitude of B, TZ is the curvature 
radius of the projected magnetic field line in XZ-plane, s measures the distance along the 
Z-axis (with the QT point at s = 0). The polar angle and azimuthal angle of B are given 
by 



The minimum 6b occurs at the QT point 

sin^t — sin 9b \s=o— e- (3.19) 
We also have 

Thus the adiabaticity parameter along the ray is 

Tad = 1.0 X 10'E,Bl,e'n, (l + , (3.21) 
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Quasi-tangential point 
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Figure 1. A sketch of the polarization mode evolution across the quasi-tangential point in the adiabatic (upper panels) and 
non-adiabatic (lower panels) limits. In both cases, the photon polarization vector E (double-arrowed bars) and the transverse 
magnetic component of the magnetic field B± (single-arrowed bars) are shown at five different positions along the ray (from 
left to right): before, slightly before, at, slightly after, after the QT point. In both limiting cases, there is no net change in the 
linear polarization when the photon traverses the QT region. Polarization change occurs only in the intermediate cases (Ft ~ 1; 
see text). 



where 7^l = 7^/(10 km). 

Before and after the QT point (where s = 0), when \s\ ^ eTZ, the adiabaticity parameter 
Tad increases rapidly with |s| and the mode evolution is generally adiabatic (Tad ^ !)• 
However at the QT point, Tad reaches its minimum value: 

Tt ~ 1.0 X lO^EiBl^e^TZi. (3.22) 

The value of Tt determines the mode evolution characteristics across the QT region. Figure 1 
shows the qualitative behaviors in two limiting cases: the adiabatic limit (Ft ^ 1) and the 
non-adiabatic limit (Tt ^ 1). In the adiabatic case, the photon polarization direction follows 
the variation of the transverse magnetic field B±. Since the final direction of B± is opposite 
to the initial direction, the final polarization direction is the same as the initial one. In the 
non-adiabatic case, ips changes rapidly, so that the polarization direction cannot follow 
and remains constant across the QT region. Thus, in both limiting cases, the polarization 
direction is unchanged when the photon traverses the QT point. 

Obviously, a non-trivial change of photon polarization across the QT point occurs only 
in the intermediate case. Ft ~ 1. To obtain the quantitative behaviors of the polarization 
evolution for general values of Ft, we integrate eq. (2.11) numerically. Figures 2-4 show 
three examples of single photon mode evolution in the QT region, corresponding to Ft <^ 
1, Ft ^ 1 and Ft ~ 1, respectively. For the Ft <^ 1 (non-adiabatic) case (Fig. 2), the 
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Figure 2. Evolution of photon polarization across the quasi-tangential region for the Ft <S 1 case. The parameters e, Ti. are 
defined in equation (3.17), s is the distance along the ray (s = corresponds to the QT point). Prior to the QT point, the 
photon is assumed to be in the ||-mode, with Ay = 1 and Ax = 0. The angles 9g, ipg specify the magnetic field orientation. 
Ax and Ay are the photon polarization amplitudes along the fixed X, Y axis, Q, U, V are the Stokes parameters. For this 
example, the parameters are B = 10^-^ G, E = IkeV, e = 5 X 10~', 1Z = 100km. Note that at the QT point, Og ~ 0, ips varies 
rapidly, giving rise to non-adiabatic mode evolution. 




close to and ipg varies slowly, giving rise to adiabatic mode evolution. 

photon polarizations (Stokes parameters) are constant throughout the region. For the Ft ^ 1 
(adiabatic) case (Fig. 3), the photon polarizations (Stokes parameters) change around the 
QT region, following the variation of B, but the final polarizations are very close to the 
initial values. The most interesting case occurs for Ft ~ 1 (Fig. 4). In this intermediate 
regime, partial mode conversion takes place, so that after crossing the QT point, the photon 
becomes a mixture of two modes (even when it is in a pure mode prior to QT crossing). The 
polarization state of the photon is therefore significantly changed by the QT effect. 

Based on the above results, we can use Ft ^ 1 to define the parameter regime for 
which propagation through the QT region gives rise to an appreciable change in the photon 
polarization. 
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Figure 4. Same as in Fig. 2, except for tlie Ft ~ 1 case, with e = 4 X 10"'^. Note that in this case, partial mode conversion 
occurs at the QT point, and the oscillatory behaviors in the Stokes parameters are the result of interference between the two 
photon modes. 

e ^ e,rit = 2.15 x 10~\E^Bl^ni)-^/\ (3.23) 

Note that when e <C ecrit, the QT effect is also neghgible. This equation effectively maps 
out the domain where the QT propagation effect must be carefully considered. Outside this 
domain, one can ignore the QT effect in determining the final observed polarization signals. 
To translate this effective parameter domain into the physical size of the emission region 
requires a knowledge of the global NS magnetic field structure. The smallness of ecrit for 
typical parameters (e.g., Ei ~ 1, S13 ~ 1 and TZi ~ 1) indicates that the "affected" region 
is a small fraction of the NS surface. But as far as the observed polarization signals are 
concerned, it is more relevant to compare the size of the "affected" region to the size of 
the photon emission area. We consider the special case of dipole magnetic field in the next 
section. 

4 QUASI-TANGENTIAL EFFECT IN DIPOLE MAGNETIC FIELD 

In this section we assume that the NS has a pure dipole field. We focus on X-ray emission 
from the polar cap region of the star (sections 4.2-4.3), but also consider more general 
emission regions on the NS surface (section 4.4). 

4.1 Magnetosphere Field Geometry Along the Ray 

To calculate the observed polarized radiation signals, we set up a fixed coordinate system 
X'Y'Z' with the Z'-axis along the line-of-sight (pointing from the NS toward the observer) 
and the X'-axis in the plane spanned by the Z'-axis and f2 (the spin angular velocity vector) 
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(a) 



(b) 



Figure 5. Geometrical model of X-ray emission in a dipole magnetic field. Photons are emitted from tlie polar cap region and 
propagate toward the observer in the direction fe along the Z or Z'-axis. At the time of emission, the magnetic dipole vector is 
/U; and the initial photon position is t\. Three coordinate systems are defined; (i) The fixed frame XYZ with Z || fc, /X; in the 
XZ-plane and /i; = /i(sin$ni;,0,cos5,ji); (ii) the fixed frame X'Y'Z' with Z \\ k, H in the X'Z'-plane and k x tl = sin^ V ; 
(iii) the instantaneous frame xyz, with z || fc, B in the x^-plane and kx B = — sinSgy. Note that as the photon propagates, 
the magnetic field it "sees" changes and thus the xyz frame rotates around the Z-axis. 

(see Fig. 5). The angle between fl and k is (. The magnetic dipole /j, rotates around O, and 
a is the angle between 17 and /i. 

Consider a photon emitted at time ti (corresponding to the NS rotation phase ^i) from 
the position ri = (i?*, Or^, (fr^) on the NS surface. At this emission time, the magnetic dipole 
moment is /ij. As the photon propagates along the Z-axis, its position vector changes as 

r = ri + sZ, (4.24) 

where s — c{t — ti) is the photon displacement from the emission point. ^ In the meantime, 
H rotates around ft, and changes according to 

IJ,{s) — fj, (sinCcosa + cosCsinacos^)^ +sinQ;sin^i^ + (cosCcoso; — sinCsinQ;cos^)Z ,(4.25) 
where the rotation phase ^ is (we set ^ = when fx lies in the X'Z' plane) 

*(s) = *i + 0(i-ii) = *i + s/nc, (4.26) 

with ric — c/Q the radius of the hght cyhnder. The polar angles {9^, ^p^) of fx in the X'Y'Z' 
frame are given by 



^ Equation 4.24 neglects the effect of gravitational light bending, which can be incorporated in a straightforward manner (e.g. 
Beloborodov 2002; van Adelsberg & Lai 2006). This effect amounts to shifting r-i in the direction perpendicular to the Z-axis, 
and does not appreciably change our result. 
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. . . > . T sin a sin , ^ 

cos^u = cosCcoso; — smCsmcKCOS w, tan(/7„ = . (4.27) 

sin C cos a + cos C sin a cos V 

(Similar expressions hold for the polar angles of /^j, Q^.^ and with ^ replaced by ^'j.) 
The changing magnetic field as "seen" by the photon is obtained from 

B{s) = -V(M • r/r^) = + ^(/x • r). (4.28) 

When discussing the polarization result of a given rotation phase (at emission) ^j, it is 
convenient to introduce another fixed coordinate system XYZ (see Fig. 5), such that Z \\ k 
and /Lt; in the XZ-plane with /U; = //(sin , 0, cos^^J. In the XYZ frame, (p/^. is given by 

fx{s) = n sin 9 1^ cos{ip ^ - ip X + sin 9^ sin^Lfi^-cp i^,)Y + cos 9 . (4.29) 

The magnetic field (4.28) as "seen" by the photon is inchned at an angle 9b with respect to 
the line-of-sight, and makes an azimuthal angle </7b in the Xy-plane. The angle 9b, <fB can 
be obtained from equation 4.28 via: 

cos9b{s)^-^, tan(^B(s) = — ^. (4.30) 
4.2 Polcirization Map of the Polcir Cap Emission 

The observed polarized radiation is the incoherent sum of photons from the emission region 
on the NS surface. For each emission point, Vi — (i?*, fn), we integrate the mode evolu- 
tion equation (2.11) along the photon ray from ri to a large radius, beyond the polarization 
limiting radius, to determine the final polarization state of the photon. The polarization 
limiting radius, Tpi, is where the two photon modes start recoupling to each other, and is de- 
termined by the condition {uj/c)An = 2\dipB/ds\. At large distance {r ^ R^,), the magnetic 
field is simply B ~ {2jjizZ — — //y^j /r^, and r^\ is given by 

'f^c.7Q{E,Bl,P,f\ (4.31) 

where B^i^ is the polar magnetic field at the stellar surface in units of 10^^ G (see van 
Adelsberg & Lai 2006 for a more detailed expression). Beyond rpi, the photon polarization 
state is frozen. As mentioned in section 1, the calculations presented in van Adelsberg & Lai 
(2006) and Lai & Ho (2003b) did not consider the possibility that the photon polarization 
may change appreciably when crossing the QT point, which typically lies at a much smaller 
radius than rpi. 

In general, the radiation emerging from the NS atmosphere at ri includes both the ||- 
mode and the ±-mode, with the intensities /ii, I±_ depending on the field strength, photon 
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Figure 6. Two-dimensional polarization map of emission from the polar cap region of a NS. The QT propagation effect 
changes the photon polarization in two narrow bands (confirmed to —Wt/2 <Y < Wt/2; see Fig. 7). The input parameters 
are: B« = 10^* G, E = IkeV and Of,. = 5° (We choose P = 5s, a = 30°, C = 35° and "I-i = 180°, although the map depends 
very weakly on these parameters as long as is the same). 



energy and emission angle (see Lai & Ho 2003b and van Adelsberg & Lai 2006). In the 
absence of the QT effect, the radiation at r > Vpi will consist of approximately the same /|| 
and I±, with a small mixture of circular polarization generated around Vpi (see van Adelsberg 
& Lai 2006). This simple result should be modified if there is a significant polarization change 
when the photon crosses the QT region. Since we are interested in understanding the QT 
effect, in the following we will assume that at the emission point the radiation is in the 
||-mode, with I\\ = 1 and I± = 0. 

Figure 6 gives an example of the 2-dimensional polarization map of the final Stokes 
parameters of photons emitted from the NS polar cap region. This map is produced for a 
specific set of parameters: surface magnetic field B^, photon energy E, and the angle between 
/Xj and k, 9^.^. [For a rotating NS, the map also depends on the angles a, the rotation period 
P and emission phase \l'i, but the dependence is mainly through 6'^. — see eq. (4.27). see 
below.] Without the QT effect, the final polarization from the polar cap region will be Q = 1. 
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Figure 7. One-dimensional polcirization profiles for emission from the polar cap region of a NS. These profiles are produced by 
fixing 6r^ while varying (pr^ of the emission point. The coordinates XY axe the projection of the emission point on the XY 
plane (see Fig. 5). The three columns (from left to right) correspond to a fixed 6^ = 4.8°, 5.0°, 5.2° (or X/R, ~ 0.837, 0.872, 
0.906), respectively. In each column, the final Stokes parameters Q, U, V axe plotted as a function of Y, as well as Ft, the 
adiabaticity parameter at the QT point. The polarizations of photons emitted from the region —Wt/2 <Y < VFt/2 (between 
the two vertical dashed lines, defined by Ft = 3) are modified by the QT effect. The input parameters for these profiles are: 
S* = 10" G, E = IkeV and e^- = 5.0° (We choose P = 5s, a = 30°, C = 35°, *i = 180°, although the profiles depend very 
weakly on these parameters as long as is the same). 



We see from Fig. 6 that the QT effect gives rise to fine features (with length scale much less 
than the polar cap size) in the polarization map. In particular, there are two narrow bands 
(which are parallel to the projected magnetic dipole axis) in the emission region, where the 
linear polarization Q is significantly reduced because the original || mode is converted across 
the QT region into a mixture of ||-mode and _L-mode. In these two narrow belts, non-zero 
pattens of Stokes U and V are produced. However, these U, V patterns arc anti-symmetric 
with respect to the Y = line (i.e., the k-fj,^ plane), so that the observed values of U and 
V almost equal zero when photons from the whole emission region are added. 

In Fig. 7 we show three sections of Fig. 6 along the Y-axis for a fixed value of X 
[X/R^ = 0.0837 or 0r, = 4.8° for panel (a), X/R^ = 0.0872 or 0r, = 9^. = 5" for panel (b), 
and X/R^ — 0.0906 or 9r^ — 5.2° for panel (c)]. We also plot the adiabaticity parameter 
Ft at the QT point for the photon rays from each emission point. The Ft profile is useful 
for understanding the results for the Stokes parameters. As the analysis in section 3 shows, 
photon propagation through the QT region changes the photon polarization significantly 
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Figure 8. Same as Fig. 7, except for = 10^* G (the left panels) and B* = 4 x lO'^^ G (the right panels), both for P = Is. 

only when Ft ~ 1. Prom Fig. 7, we see that the photon Stokes parameters are modified by 
the QT effect only in two narrow bands, confined in the region —Wt/2 < Y < Wt/2, with 
the boundaries defined approximately by Ft = 3 (the two vertical dotted lines). Outside 
this region, Ft ^ 1, the mode evolution is adiabatic around the QT point, so that the 
photon polarization state follows the magnetic field direction and the Stokes parameters are 
unchanged across the QT region. In the middle of the band (i.e., around y = 0), Ft <S 1, 
mode evolution is non-adiabatic, and the photon polarization state also doesn't change across 
the QT point (see Fig. 1). 

Figure 8 depicts other examples of the ID polarization map for different values of stellar 
magnetic field (10^"^ and 4 x 10^^ G) and spin period P. As in Fig. 7, we see that the final 
photon polarizations are determined by the value of Ft, the QT effective region is confined 
to —Wt/2 < Y < Wt/2 defined by Ft < 3 (the region between the two vertical dotted fines). 

A careful examination of Figs. 7-8 shows that due to the NS rotation, the Ft profile 
and Stokes profiles shift in Y by the amount AY ~ (i?^, sin Q;/ric)st, where St is the distance 
of the QT point from the emission point. Since typically St is less than a few NS radii, this 
shift AY is negligible. According to eq. (3.22), Ft is proportional to EB^, thus for smaller 
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Figure 9. The width of the QT effective region — VKt/2 <Y < Wt/2, defined by Ft < 3, as a function of the surface magnetic 
field strength and photon energy. Here B«i4 = B«/(10^* G), Ei = B/(lkeV). The different lines correspond to different 0^, : 
e^. = 5° (soUd Une), 6^, = 15° (dashed Une), 0^^ = 45° (dot-dashed line). 



-B* (and smaller E) the effective width Wt is larger, as seen in Fig. 8. In Fig. 9 we show how 
Wt changes with varying EB^. Our numerical result for Wt can be fitted by 

1^ ^ 2.7 X 10-^ iBl,,E,)-'/'fi9,,). (4.32) 

Here f{0^^) is a dimensionless function of [which in turn depends on a, C, and ^'i, see 
eq.(4.27)]: f{9^^ = 5°) = 1 and varies from 1 to 1.7 for different values of (see Fig. 10). 
The scaling relation in eq. (4.32) is derived in the appendix. 
Note that the polar cap width (diameter) is given by 



^^cap ~ / R* 



R* \ c/Q, 
where P5 = P/(5s). Then we have 



1.3 X lO-^Ps"^/^ (4.33) 



^ = 2.1{Bl,,E,P^-"')-''^f{e,.:). (4.34) 



W^cap 



Equation (4.34) implies that the size of the effective QT region (where the QT effect changes 
the photon polarization) can be comparable to the polar cap size for some parameters (e.g., 
low photon energy and low magnetic field strength). 
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Figure 10. The dimensionless function f{6^n), as defined in eq. (4.32). 



4.3 Quasi- Tangential Effect on the Observed Polarization 

As discussed above, when linearly polarized radiation with Q — I — 1 traverses the QT 
region with Ft < 3, its Stokes parameters will be changed (so that Q will become less than 
unity, and V will be nonzero; see Figs. 7-8). However, when adding up radiation from 
a finite emission region, we find / dYU ^ and / dYV ~ 0. Thus the net effect of the 
QT propagation is to reduce the degree of the linear polarization of the photon. In general, 
if Fq is the flux of linearly polarized radiation prior to passing the QT region, then after 
traversing the QT region, the linearly polarized radiation flux Fq can be obtained by 



where Wem is the width of the emission region. As seen in section 4.2 (see Figs. 7-8), the 
Q{Y) proflles for different parameters are similar, so that Fq/ Fq depends only on Wt/Wem- 
Our numerical result for Fq/ Fq as a function of Wt/Wcm is shown in Fig. 11. For a given 
set of parameters {B^, E, ^^.), we can use eq. (4.32) and Fig. 10 to obtain Wt, and then 
use Fig. 11 to read off the ratio Fq/Fq. If ^ Wem, we flnd that Fq/Fq is approximately 
given by 




(4.35) 



Fq/Fq ~ 1 



(for m/Wem < 1). 



(4.36) 
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Figure 11. The effect of QT propagation on ttie observed linear polarization: Fq (Fq) is the linearly polarized radiation flux 
before (after) traversing the QT region; see eq. (4.35). Here Wt is the width of the effective QT region (where Ft < 3) and 
Won is the width of the emission region. For Wt/Wcm ^ 1, Fq/Fq can be approximated by eq. (4.36), which is shown as the 
dashed line. 



Note that for <^ W^em? the effective QT region is much smaller compared to the emission 
region so that Fq/Fq ~ 1; for Wem <^ Wt, the photon mode evolution across the QT 
region is non-adiabatic so that Fq/Fq is also close to unity. Thus in both Wt/W <^ 1 and 
Wt/Wcm ^ 1 limits, Fq/Fq = 1, i.e., the linear polarization is unchanged by the QT effect. 
The minimum value of Fq/Fq occurs at PVt/W^m — 1-8. 

Suppose the emission is coming from polar cap region with width Wcm = W^cap, then 
according to eq. (4.34), Wt/Wcap is proportional to (Bli^Ei)'^^^ P^^'^ . Thus, for given -B=k13, 
P5 and 6^., we can translate Fig. 11 into a spectrum for Fq/Fq — This spectrum is shown 
in Fig. 12. We see that at sufficiently low energies and high energies, Fq/Fq ^ 1, and 
the minimum reduction of the degree of linear polarization due to the QT effect occurs at 
Wt/Wca.p — 1.8, corresponding to the photon energy 



Now we consider the (energy-dependent) polarization light curve produced by the polar 
cap emission of a rotating NS. As the NS rotates, the k-fj,^ angle 9^^^ changes with the 
rotation phase [see eq. (4.27)]. According to eq. (4.34) and Figs. 10, 11, different 9fj_. will 



^min ^ 1.6B;,lPrfiO^,)keV. 



(4.37) 
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Figure 12. The reduction factor Fq/Fq of the linear polarization as a function of the photon energy (and other parameters). 
Here -B*i3 = -B,/(10^^ G), -Ei = -E/(lkeV), P5 = P/(5s), and we have assumed V14m = VKcap. The angle between A; and /X; is 
chosen to be = 5°. 

give different Fq/Fq , which means Fq/Fq will evolve with the rotation phase \E'i. Moreover, 
Wt/Wcap and Fq/Fq depend on the photon energy. Figs. 13 and 14 present two examples 
of the phase evolution of Wt/Wcap and Fq/Fq, for = 10^'^ G and 10^^ G, respectively. 
As noted before, Fq/Fq reaches a minimum at Wt/Wca.p = 1-8 (see Fig. 11). In the case 
depicted in Fig. 13, Wt/Wcap > 1-8 for E = IkeV and 0.5 keV at all phases, thus Fq/Fq is 
larger for lower photon energies (since Wt/Wca.p increases with E). In Fig. 14, W^/Wcap < 1-8 
for all the three photon energies (i? = 0.5, 1, 5keV), so that Fq/Fq is larger for higher 
energies. 

To produce the observed polarized radiation fluxes, using the results present above, we 
will need input Fq. In general, the emergent radiation from the NS atmosphere (before 
passing through the QT region) is linearly polarized, with Fq dependent on E and the 
rotation phase (and a, Q. In the XY Z frame, the radiation has polarized fluxes Fq ^ 
(which can be either positive or negative), Fu = and Fy = 0. Using Fq from an atmosphere 
model and our result for Fq/Fq, we can compute Fq, the polarized flux after passing through 
the QT region. Then the observed polarized radiation fluxes in the fixed X'Y'Z' frame are 
given by 
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Figure 13. The phase evolution of the modification of Unear polarization by the QT effect, F'q/ Fq (upper panel) and the 
effective QT region width relative to the polar cap size, W/Wcs.p (lower panel). The different curves are for different photon 
energies: E =0.5 keV (dashed lines), IkeV (solid lines) and 5keV (dotted line). The input model parameters are: _B* = 10^'^ G, 
P = 5s, a = 30° and C = 35°. The emission region size is assumed to be Wcap- 
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Figure 14. Same as Fig. 13, except for = 10^'* G. 
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F'q = cos 2(^s (rpi) , F{j = Fq sin 2(pB (rpi) , 



(4.38) 



where </?B(rpi) is the azimuthal angle of the magnetic field, (ps, evaluated at the polarization 
limiting radius r = Tpi. For rpi <^ r\c = c/fl, or for spin frequency <S 70 Hz, (/^^(rpi) ~ 
TT + ^ TT + (see Lai & Ho 2003b; van Adelsberg & Lai 2006), thus Fq = Fq cos 299^. 
and FIj = Fq sin 2(pf,.. 

4.4 Quasi- Tangential Effect for Different Emission Regions 

In the previous subsections we have focused on emission from the polar cap region of the 
NS. In reahty, the "hot spot" on the NS surface may be larger. The QT effect also exists 
outside the polar cap region. 

As discussed above, the QT effective region is confined to Ft < 3. For a given dipole field 
geometry, we can calculate Ft for photon rays emerging from different points on the stellar 
surface. Figs. 15 and 16 give two examples of the effective QT region as defined by Ft < 3, 
for the fe-/Ui angle 9^. — 15° and 5", respectively. For each emission point, if we ignore the QT 
effect, the final photon polarization is dominated by the wave mode coupling effect and the 
polarization position angle (PA) is approximately determined by (pB at polarization limiting 
radius rpi (again, assuming at emission the photon is in the ||-mode). Generally we have 
R* ^ Tpi <^ ric, so (pB{rp\) ~ (p^. +7r. Therefore, the final polarization direction of the X-ray 
photon from any point of the star surface are always parallel to the plane. However, in 
the QT effective region, the final PA may be modified when the photon propagates through 
QT region (see the lower panel of Fig. 15). Note that the final polarization angles are 
always symmetric with respect to the plane. Thus, as discussed before, when adding 
up radiation from a finite emission region, the net effect of QT propagation is to reduce the 
degree of linear polarization without changing the polarization angle. 

Special caution must be taken when is small. We see from Fig. 16 that for — 5", 
the effective QT region may cover a significant part of the stellar surface (In Fig. 16, the 
hatched region only shows the effective QT region that satisfies < 20".) This can be 
understood as follows. For i?* <^ r <S r\c the transverse parts of the dipole magnetic field 
are [see eq. (4.28)] 



with r — ri + sZ (neglecting light bending) and s ^ ri, and Xi, Yi are the X, Y components 




(4.39) 
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of the emission position Vi. For 6*^. ^ 3i?*/r, we have Bx — — (/x/r^) sin6'^., By — 0, thus 
outside the QT region, the final polarization is aligned with the k — /j,^ plane. However, for 
Ofj,. <, 3i?*/r, we no longer can neglect By relative to Bx- Thus for sufficiently small 
polarization alignment will not be achieved for most emission points. Figure 17 shows the 
ID profiles of the final Q, U, V and Ft for a fixed and varying (p^, all with 9^. — 5°. 
In the left and middle panels where 6',.. = 10° and 20°, the profiles are similar to the cases 
examined before (see Figs. 7 and 8). For the Oj.^ = 30° case (the right panels of Fig. 17), 
Ft is always less that 3. The reason is that when 9^ is large, the QT point hes far away 
from the star rt ^ R* (and it is even possible that rt > Vpi for sufficiently large 9^)- In this 
case (small 9i^. and emission from the region far away from the magnetic pole), the simple 
prescription we have presented in section 4.2 - 4.3 to account for the QT effect cannot be 
used, and numerical ray integration from each emission point is necessary. 



5 DISCUSSION 

We have studied the evolution of photon polarization in a neutron star magnetosphere 
whose dielectric property is dominated by vacuum birefringence. We have focused on X-rays 
because of the potential of using X-ray polarimetry to constrain neutron star magnetic fields 
and to probe strong-field QED (see Heyl & Shaviv 2002; Heyl et al. 2003; Lai & Ho 2003b). 

5.1 X-ray Polcirization Signals without QT Effect 

If one neglects the QT propagation effect studied in this paper, then it is straightforward to 
obtain the obserevd polarized X-ray fluxes (Stokes parameters) from the fluxes at the emis- 
sion region, at least approximately, without integrating the polarization evolution equations 
in the magnetosphere (Lai & Ho 2003b; van Adelsberg & Lai 2006). For a given (small) 
emission region of projected area^ ^A±, one need to know the intensities of the two photon 
modes at emission, I± and In the case of thermal emission, these can be obtained directly 
from atmosphere/ surface models. As the radiation propagates through the magnetosphere, 
the photon mode evolves adiabatically, following the variation of the magnetic field, until the 
polarization limiting radius Tpi, at which point the polarization is frozen. Thus the polarized 



^ This is the area perpendicular to the ray at the emission point — General Relativistic light bending effect can be easily 
included in this. 
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Figure 15. The quasi-tangential effect for X-rays emitted from tfie whole surface of the NS with the k-n^ angle 9^^ = 15°. 
Here we use the XY Z frame with fc || Z, /j.; in XZ-plane and the circle corresponds to the project star surface. In the upper 
panel, the hatched region is the effective QT region, defined by Ft < 3. The bars show the polarization direction. The cross 
marks the projected dipole axis ^x^ on the stellar surface. The lower panel shows the enlarged effective QT region. The other 
input parameters are: -B» = 10^* G, E = 1 keV. 

radiation flux beyond r^x is Fq = (J|| — Ij_)AAj_/D, and Fu ~ Fy ~ 0^, where D is the 
distance of the source, Fq and Fu are defined in the coordinate system such that the stellar 
magnetic field at Vpi lies in the XZ plane (with the Z-axis pointing towards the observer). 
Since Vpi is much larger than the stellar radius, the magnetic fields as "seen" by different 
photon rays are aligned and are determined by the dipole component of the stellar field, one 
can simply add up contributions from different surface emission areas to Fq to obtain the 
observed polarization fluxes. 

Note that the description of the polarization evolution in the last paragraph is valid 

^ Note that Fy is not exactly zero because of the neutron star rotation and because mode recoupling does not occur instantly 
at rpi; see van Adelsberg & Lai (2006). 
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Figure 17. Similar to Fig. 7, but for emission from larger regions of the NS. Here we use the varying as the X-axis. The 
three columns (from left to right) correspond to a fixed 9^. = 10°, 20°, 30° (or X/R, = 0.174, 0.342, 0.5 at Y = 0). The 
input parameters are = 10^* G, E = IkeV, 9^^. = 5°. 
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regardless of the possible complexity of the magnetic field near the stellar surface. This 
opens up the possibility of constraining the surface magnetic field of the neutron star using 
X-ray polarimetry. For example, the polarization light curve (particularly the dependence on 
the rotation phase) depends only on the dipole component of the magnetic field, while the 
intensity lightcurve of the same source depends on the surface magnetic field. On the other 
hand, the linear polarization spectrum (i.e., its dependence on the photon energy) depends 
on the magnetic field at the emission region (Lai & Ho 2003b; van Adelsberg 2006); thus it 
is possible that a NS with a weak dipole field (^ 7 x 10^^ G) may exhibit X-ray polarization 
spectrum characteristic of a ^ 7 x 10^^ G NS (see section 1). 



5.2 Effect of QT Propagation 

The QT propagation effect studied in this paper complicates the above picture somewhat. 
When a photon passes through the QT region (which typically lies within a few stellar radii) , 
its polarization modes can be temporarily recoupled, and this can give rise to partial mode 
conversion. Thus after passing through the QT point, the mode intensities change to /|| ^ I\\ 
and Ia_. The observed polarization flux is then Fq = (/|| — I±)^A/ D ^ Fq. 

In the most general situations, to account for the QT propagation effect, it is necessary 
to integrate the polarization evolution equations (see section 2.2) in order to obtain the 
observed radiation Stokes parameters. However, we show in this paper that for generic 
near-surface magnetic flelds, the effective region where the QT effect leads to signiflcant 
polarization changes covers only a small area of the neutron star surface [see section 3, 
particularly eq. (3.22)]. For a given emission model (and the size of the emission region) and 
magnetic field structure, one can use our result in section 3 to evaluate the importance of 
the QT effect. In the case of surface emission from around the polar cap region of a dipole 
magnetic field, we have quantified the effect of QT propagation in detail (section 4) and 
provided a simple, easy-to-use prescription to account for the QT effect in determining the 
observed poUarizatin fluxes. Our key results are presented in sections 4.2-4.3, particularly 
Fig. 11 and equations (4.32) (together with Fig. 10) and (4.35) - (4.36). As discussed in 
section 4.3, the net effect of QT propagation is to reduce the degree of linear polarization, so 
that Fq/Fq < 1, with the reduction factor depending on the photon energy, magnetic field 
strength, geometric angles, rotation phase and the emission area. The largest reduction is 
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about a factor of two, and occurs for a particular emission size (VFt/W^cm — 1-8; see Fig. 11). 
Obviously, for emission from a large area of the stellar surface, the QT effect is negligible. 

Overall, the QT effect will have a small or at most modest effect on the observed X-ray 
polarization signals from magnetized neutron stars. Thus the the description of X-ray polar- 
izaton signals given in section 5.1 remains largely valid. In most situations (the exceptions 
are discussed in section 4.4), the QT effect can be accounted for in a straightforward manner 
using the result presented in our paper. 
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APPENDIX A: DERIVATION OF EQUATION (4.31) 

Consider the emission from the polar cap region of a non- rotating NS. In the XY Z-lrome. 
(see Fig. 5, k || Z, fx^ in XZ-plane), the initial magnetic momentum is 

/Ui = //(sin6l^.,0,cos6l^;), (Al) 

and the photon position is 

r — {R^ sin 9j.. , Y, cos + s) , (A2) 

where s is the displacement from the emission point, and Y <^ i?^. We want to evaluate Ft 
at the QT point s — St){Y) as a function of Y. The magnetic field is given by eq. (4.28), 
and has components: 

„ sin du- 3i?* sin 9n. , ^ ^ , 

Bx = ^5 '' {R* + s cose,, I 

?>Y 

By = —{R^ + s cose,.), (A3) 

(here we set ^ = 1). For a given F, The QT point (s — Sx{Y)) is defined by 

dBi 

ds ^~ 

We can expand Bx and By near s — s^iY) using Taylor expansion: 



0. (A4) 



Bx = Bxt + FiAs 

By = {H + F2As)Y (A5) 
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with As = s — Sf Note that Bxt, H, Fi, F2 all depend on St{Y). Substitute eq. (A5) into 
eq. (A4), we have 

B., = -r^:^. (A6) 

From eq. (A6) we see that 

As,(Y)^st{0) + O(Y^). (A7) 
Thus 

F^ = Fio + 0{Y'), F, = F^o + 0{Y'), H^Ho + 0{Y'), Bxt = -Y'^^ + 0{Y') (A8) 

where Fio, F20, Hq arc constants. . The azimuthal angle of B is determined by tan(^B = 
By/Bx- Therefore ip'^ at s = St is 

/ , 2 fBvV , YHFi + Y^HFi/Fi , 

Is^s. = COS ^— j |_ = oc l/Y (A9) 

We also have 

sm9B\s=s,= ^<xY (AlO) 

Substitute eq. (A9) and (AlO) into eq. (2.13), the adiabaticity parameter at QT point is 
then given by 

Ft oc BlEY\ (All) 
The QT effective width Wt, defined by Ft = 3, is then Wt oc {B^E)''^/^. 



